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Abstract 

A Kondo-lattice theory is applied to the t-J model on a quasi-two dimensional lattice. The 
Kondo temperature Tk or /cbJk is defined as an energy scale of local quantum spin fluctuations, 
with k-Q the Boltzmann constant. The bandwidth W* of quasiparticles, which is approximately 
given by W* ~ ^k^T^i, is renormalized by the Fock-type term of the superexchange interaction. 
When the lifetime width 7 of quasiparticles is so small and temperature T is so low that 7/VF* ^ 1 
and k-^T /W* ^ 1, the renormalization is large. Even in the limit of the half filling, for example, the 
bandwidth W* is nonzero and of the order of | J|, with J the superexchange interaction constant 
between nearest neighbors. The Kondo temperature /cbTk, which also gives a measure of the 
strength of the quenching of magnetic moments by local quantum spin fluctuations, increases 
away from the half flling; it is smaller for larger 7. Therefore, local-moment magnetism, which is 
characterized by Tn ^ Tk, with Tn the Neel temperature, appears for almost half flllings, and 
itinerant-electron magnetism, which is characterized by Tn ^ Tk , appears for fillings away from the 
half filling; Tn is higher for larger 7. The antiferromagnetic region as a function of electron fillings 
is wider for larger 7; this result implies that it is wider in more disordered systems than it is in less 
disordered ones. The difference or asymmetry of disorder between electron-doped and hole-doped 
cuprate oxide superconductors must be, at least partly, responsible for that of antiferromagnetic 
phases between them. 
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I. INTRODUCTION 



The discovery^ in 1986 of high critical-temperature (high-Tc) superconductivity in cuprate 
oxides has revived intensive and extensive studies on strong electron correlations because 
it occurs in the vicinity of the Mott-Hubbard metal-insulator (M-I) transition or crossover. 
When no electrons or holes are doped, cuprates are insulators; they exhibit antiferromag- 
netism at low temperatures. When enough electrons or holes are doped, they become metals; 
they exhibit superconductivity at low temperatures. The M-I transition or crossover and 
the transition between antiferromagnetic and paramagnetic phases as a function of doping 
concentrations are asymmetric between electron-doped and hole-doped cuprates^; the insu- 
lating phase is much wider in electron-doped ones than it is in hole-doped ones, and the 
Neel temperature is relatively higher in electron-doped ones than it is in hole-doped ones. It 
is plausible that a relevant theory can explain not only high-Tc superconductivity but also 
the asymmetry of the M-I transition or crossover and that of the Neel temperature. 

In 1963, three distinguished theories on electron correlations in a single-band model, 
which is now called the Hubbard model, were published by Kanamori,'^ Hubbard,— and 
Gutzwiller.— Two theories among them are related with the M-I transition or crossover. 
According to Hubbard's theory— the band splits into two subbands called the lower and 
upper Hubbard bands. According to Gutzwiller's theory,^ with the help of the Fermi-liquid 
theoryi^ a narrow band of quasiparticle appears on the chemical potential; we call them 
Gutzwiller's band and Gutzwiller's quasiparticles. The combination of the two theories 
imphes that the density of states must be of a three-peak structure, Gutzwiller's band 
between the lower and upper Hubbard bands. This speculation was confirmed in a previous 
paper.— The Mott-Hubbard splitting occurs in both metallic and insulating phases as long 
as the onsite repulsion U is large enough, and Gutzwiller's band is responsible for metallic 
behaviors. Then, we can argue that a metal in the vicinity of the M-I transition can become 
an insulator only when a gap opens in Gutzwiller's band or that it can behave as a bad 
metal or an insulator when lifetime widths of Gutzwiller's quasiparticles are larger than 
their bandwidth. 

Not only Hubbard'a^ and Gutzwiller'^ theories but also the previous theory^ are within 
the single-site approximation (SSA). This fact implies that local electron correlations must be 
responsible for the three-peak structure. Local electron correlations are rigorously considered 
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in SSA that includes all the single-site terms; such an SSA is rigorous for paramagnetic 
phases with no order parameter in infinite dimensions.— It is reduced to solving the Anderson 
model /°i^^i^^ which is one of the simplest effective Hamiltonians for the Kondo problem. 
This approximation is often called the dynamical mean- field theory (DMFT).J^ The Kondo 
problem has been almost settled or almost completely clarifie d^"^i'^^i'^^i^'''i^^i^^ . The Kondo 
temperature Tk or A;b^k is defined as an energy scale of local quantum spin fiuctuations, 
with k-B the Boltzmann constant. Gutzwiller's band between the lower and upper Hubbard 
bands corresponds to the so called Kondo peak between two sub-peaks; its bandwidth W* 
is about W* ~ AksTK- 

One of the most essential physics involved in the Kondo problem is that a magnetic mo- 
ment is quenched by local quantum spin fiuctuations so that the ground state is a singlelji^ or 
a normal Fermi liquid . A strongly correlated electron liquid on a lattice is a normal Fermi 
liquid at T ^ Tk, where magnetic moments are quantum- mechanically quenched, while it is 
a nondegenerate Fermi liquid at T ^ Tk, where magnetic moments are thermally quenched. 
Local-moment magnetism occurs at T ^ Tk, and itinerant-electron one occurs at T ^ Tk; 
superconductivity can only occur at T <^ Tk. The M-I transition or crossover is related 
with the crossover between local-moment magnetism and itinerant-electron magnetism. 

The coherence of quasiparticles is destroyed by not only thermal fiuctuations but also 
disorder. When ^ W* or •y ^ W*, with 7 the quasiparticle lifetime width, quasi- 
particles are never well-defined; the liquid behaves as a bad metal or an insulator. When 
kBT ^ W* and 7 ^ W*, on the other hand, they are well-defined; the hquid behaves as a 
metal. Disorder can also play a significant role in the M-I transition or crossover. 

Kondo-lattice theoriest22iSii2^ have been developed so far mainly in order to elucidate 
strong electron correlations in typical or canonical Kondo lattices of lanthanide and actinide 
alloys, which are also called dense Kondo systems or heavy-fermion systems. In these canon- 
ical Kondo lattices, almost localized 4/ or 5/ electrons and itinerant conduction electrons 
coexist and conduction electrons can play an important role not only in the quenching of 
magnetic moments of almost localized electrons but also in the appearance of an intersite ex- 
change interaction mediated by conduction electrons such as the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) exchange interaction^^^ and the opening of the so called hybridization gaps or 
pseudogaps. In this paper, we are interested in electron correlations in non-canonical Kondo 
lattices such as transition-metal alloys in the vicinity of the M-I transition or crossover. 
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where the coexistence of almost locahzed electrons and conduction electrons is absent or it 
plays no essential role even if it is present. 

A Kondo-lattice theory developed in this paper is a perturbation theory starting from 
a starting or an unperturbed state constructed in a non-perturbative method, that is, 
SSA or DMFT; intersite terms are perturbatively considered in terms of intersite ex- 
change interactions. This perturbation is nothing but 1/d expansion, with d being spa- 
tial dimensionality.— liSiii It has already been applied to typical issues on electron corre- 
lations in non-canonical Kondo lattices such as the Mott-Hubbard M-I transition,— the 
Curie- Weiss law of itinerant-electron magnets,— itinerant-electron antiferromagnetisnt^ 
and ferromagnetism,— and related issues to high-Tc superconductivity,— and so on. Early 
paper o^°i^^ on a mechanism of d'-f-wave high-Tc superconductivity can also be regarded within 
the framework of the Kondo-lattice theory. The interplay between correlations and disorder 
in a paramagnetic phase is investigated for the Hubbard model with the just half filling.— 
One of the purposes of this paper is to study the interplay between correlations and dis- 
order in the crossover between local-moment magnetism and itinerant-electron magnetism 
as a function of electron fillings in the t-J model. This paper is organized as follows: The 
Kondo-lattice theory is reviewed in § [H] A paramagnetic unperturbed state is constructed 
in ^ IIIII Antiferromagnetic instability of the unperturbed state is studied in S IIVI Discussion 
is given in § Conclusion is given in § IVII The selfenergy of quasiparticles in disordered 
Kondo lattices is studied in Appendix. 

II. KONDO-LATTICE THEORY 
A. Renormalized SSA 

We consider the t-J mode l??i?^ or the t-t'-J-Uoo model with Uao/\t\ +00 on a quasi-two 
dimensional lattice composed of simple square lattices: 

n = -tY^ alaj^ -t'Y^ alaj^ - ^ J ^(5, ■ S^) + U^Yl rii^nn, (2.1) 

{ij)o- (ij)'a (ij) i 

with t the transfer integral between nearest neighbors (ij), t' between next-nearest neighbors 
{ij)' , Si = J2ai3 I {'^x^ ^ ^y'^ ^ ^z^) '^la'^i^^ with (7^;, (Jy, and (Jz being the Pauh matrices, and 
nia = al^aicr- Infinitely large onsite repulsion Uoo is introduced to exclude doubly occupied 
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sites. The t-J model ()2.1|) itself can only treat less-than-half fillings. When the hole picture 
is taken, it can also treat more-than-half fillings. We consider two models: a symmetric one 
with t' = 0, where physical properties are symmetric between less-than-half and more-than- 
half fillings, and an asymmetric one with t' 7^ 0, whose precise definition is made in § 11111 
Effects of disorder and weak three dimensionality are phenomenologically considered in this 
paper. We assume J /\t\ = —0.3; t = (0.3-0.5) eV and J = —(0.10-0.15) eV for cuprate oxide 
superconductors . 

We follow the previous paper— to treat Uoo- The selfenergy is divided into the single-site 
one T,„{ien) and the multisite one ASo-(ie„, k). The single-site one T,„{ien) is given by that 
of the Anderson model with the same Uoo as that of the t-J model (j2.H) : other parameters 
should be selfconsistently determined to satisfyiSJi 

G'.(i5„) = ^$^a(ie„,fc), (2.2) 

k 

with the number of sites, and G'o-(i^n) and 

G^iiEn, k) = ^ (2.3) 

ien+Ai--E(fc)-S^(i£:„)-AS^(i£:„, k) 

the Green functions of the Anderson and t-J models, respectively. Here, n is the chemical 

potential and E{k) = —2triis{k) — 2t'r]2s{k), with 

Vis{k) = cos{kxa) + cos{kya), (2.4) 



V2s{k) = 2cos{kxa) cosikyO). (2.5) 

with a the lattice constant of the simple square lattices. 

Note that the mapping condition ()2.2|) with eq. ()2.3|) depends on the multisite selfen- 
ergy, which should be perturbatively and selfconsistently calculated. There can be various 
approximations of singles-site selfenergies or dynamical mean fields depending on it. As 
is shown in § 11111 the multisite selfenergy includes an energy independent part, which is 
denoted by AEo-(fc), because the superexchange interaction constant does not depend on 
energies in the t-J model. In this paper, we include only /\'Lcr{k) in eq. ()2.3|1 : the Green 
function that includes ASo-(fc) instead of ASa.(ie„, k) is denoted by G^a\\Sn, k). 

The Kondo temperature Tk is defined by 

A;bTk = [1/Xs(0)]2.^0' (2-6) 
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where Xs(ii^i) is the spin susceptibihty of the Anderson model, which does not include the 
conventional factor (l/4)(7^/i|, with g and /ie being the g factor and the Bohr magneton, 
respectively. The single-site selfenergy T,„{ien) is expanded for < 2A;b^k as Sa-(i£^n) = 
S(0) + (1 — (p^yiSn + Ylia'^^ ^ (paa')^fJ'a' + ' ' ' 5 wlth (f)^ = (j)^^ and /S.^„ infiniteslmally 
small spin-dependent chemical potential shifts. The Wilson ratio is defined by Ws = 4>s/4'-y, 
with 0s = 0cr<T — (pa-a- FoY almost half fillings, charge fluctuations are suppressed so that 
4>c = 4>cTa+4>cT-a- ^ 1- For such fiUings, 0^ 1 so that — 2. The Green function is divided 
into coherent and incoherent parts: G^a\i£n, k) = {l/(f)^)ga\ien, fc)-|-(incoherent part), with 

k) = . ^ , (2.7) 

with fi* = {fi — Eo)/07 and ^(fc) = [E{k) + AS(fc)]/0^. We introduce a phenomenological 
lifetime width 7. Although 7 depends on energies even if it is due to disorder, as is discussed 
in Appendix El its energy dependence is ignored.— 

According to the Fermi-surface sum rule,-^ the number of electrons per site for T/Tk — >■ 
+0 and 'y/k^TK +0 is given by 

n = 2 y dep^-,o{e)f^{e - /i*) = 2 J dep^{e)fj=o{6 - /i*), (2.8) 

with 

P^(^^ = V]^T.} ,,l2^ (2-9) 



Z n 



(2.10) 



2 2nkBT^ 

with ip^z) the di-gamma function. Note that f^=o{e) = 1/ [exp(e/A;B7') + 1]. We assume 
eq. ()2.8p even for nonzero T and 7. The parameter Sq or p* can be determined from 
eq. ()2.8|) as a function of n. 



B. Intersite exchange interaction 

In Kondo lattices, local spin fluctuations at different sites interact with each other by 
an exchange interaction. Following this physical picture, we define an intersite exchange 
interaction Is{i<^i, q) by 

Xs(i^i, q) = TTT V^J-^^ (2.11) 

1 - jls{iuji, q)Xs{i(^i) 



where Xs{^<-^i,q) is the spin susceptibihty of the t-J model to be studied. Following the 
previous paper,— we obtain 

Uiui, q) = J{q) + 2f/^A7rs(iu;/, q), (2.12) 

with J{q) = 2Jriis{q) and ATTs{iuJi, q) the multi-site part of the irreducible polarization 
function in spin channels; 2U^ATTs{iuJi, q) is examined in ^ II V Bl and ^ IIV CI 

According to the Ward relation,™ the irreducible single-site three-point vertex function 
in spin channels, \s{ien,i£n + ^<^h^<-^i), is given by 

UooK{i£n, i£n + i^^i; iuji) = 24>Jxs{m), (2.13) 

for |e„| — > +0 and \u!i\ +0. We approximately use eq. ()2.13p for < 2A;b^k and 
^ 2A;b2k- A mutual interaction between quasiparticles mediated by spin fluctuations is 
given by 



203/xs(i^O ' i^s(i^/, q) = 0s^4*(i^^ Q), (2-14) 
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with 

Fsiii^h q) = Xs(it^i, q) - Xs(it^i)> (2-15) 

^T*f \ i^s(i^«,q) 

-I, {lui, q) = . (2.16) 

The single-site term is subtracted in Fs{iuJi, q) because it is considered in SSA. The expansion 
coefficient 0s appears as an effective or reducible single-site vertex function in eq. ()2.14p . 
Because the spin space is isotropic, the interaction in the transversal channels is also given by 
these equations. In the Kondo-lattice theory, intersite effects are perturbatively considered 
in terms of hiy^j^u q)i 4*(i^«) 9); or Fsiioji, q) depending on each situation. 



III. RENORMALIZATION OF GUTZWILLER'S QUASIPARTICLES 

In this section, we restrict the Hilbert space within paramagnetic states. There are two 
types of selfenergies linear in the exchange interaction: Hartree-type and Fock-type terms. 
The Hartree-type term vanishes in paramagnetic states. We consider the Fock-type term of 
the superexchange interaction: 

AS.(fc) = ^02:fc^5^J(fc-p)e-"°U(ien,p). (3.1) 
7 
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FIG. 1: Renormalized transfer integrals t* of the symmetric model: (a) kBT/\t\ = 0.02, (b) 
kBT/\t\ = 0.1, (c) kBT/\t\ = 0.2, and (d) kBT/\t\ = 0.4. In each figure, topmost solid, dashed, 
broken, chain, and chain double-dashed lines show results for 'y/\t\ = 0.01, 0.1, 0.2, 0.4, and 1, 
respectively. For comparison, l/(p.y is also shown by a bottom solid line. 



The factor 3 appears because of three spin channels. When only the coherent part is con- 
sidered, 

1 q _ 

(3.2) 



CD., 4 



with 



(3.3) 



The dispersion relation of quasiparticles is given by 



^{k) = -2t*r]Uk)-2t;r]2s{k), 



where t* should be selfconsistently determined to satisfy 



2t* = ^- -W^JE, 



(3.4) 



(3.5) 



and ^2 is simply given by = t' /(p^. 

For the symmetric model, t' = so that = 0. In order to examine how crucial role 
the shape of the Fermi surface plays in the asymmetry, we consider a phenomenological 
asymmetric model with t2/t* = —0.3. 

Expansion parameters (p^ and 0s are given by those of the mapped Anderson model. 
However, we approximately use those for the Anderson model with a constant hybridization 
energy. According to Appendix of the previous paper—. 



1/1 A (vr/2)^(l-|^|)^ 
2V|5| 'V cos2(7rV2) ' 



(3.6) 
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FIG. 2: t* of the asymmetric model for k-QT/\t\ = 0.02. See also the caption of Fig. ^ 
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FIG. 3: Single-particle properties of the unperturbed state of the symmetric model, (a) Density 
of states for quasiparticles /0^(e), (b) effective chemical potentials, fi* , as functions of carrier con- 
centrations n, and (c) n as functions of /x*. In these three figures, solid, dashed, broken, chain, and 
chain double-dashed lines show results for 7/|t*| = 10"'^, 0.1, 0.3, 1, and 3, respectively. Note that 
the solid and dashed lines are the almost same as each other, (d) Fermi surfaces for k^T /\t\ = 0, 
7 = 0, and 19 electron concentrations such as n = 0.1 x i, with 1 < i < 19 being an integer. Dashed 
lines show Fermi surfaces for n =0.5, 1.0, and 1.5. Single-particle properties of the asymmetric 
model can be found in Fig. 2 of ref. |28 . 



~ _ 1 {iT/2Y{l-\5\f 

where 5 = 1 — n is the concentration of dopants, holes [5 > 0) or electrons (5 < 0). These 
are consistent with Gutzwiller's theory.— Figures ^ and El show t* of the symmetric and 
asymmetric models as a function of 6 for various 7. It is interesting that t* is nonzero even 
for 5 if 7 and k^T are small enough. For the symmetric model (tg = 0), eq. (jH.Hj) can be 
analytically calculated for T = K, 7 = and 6 = (/i* = 0) in such a way that H = A/tt"^. 
It follows that [t*/t]g^Q - (3W^J/2nH^ = 0.18 for J/t = -0.3. If 7 or fceT is large 
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enough, on the other hand, S and t* vanish for S ^ 0. 

Figure 01 shows physical properties of the unperturbed state of the symmetric modeL 
Those of the unperturbed state of the asymmetric model can be found in Fig. 2 of ref. |28. 



In § IIVI we study antiferromagnetic instabilities of these paramagnetic unperturbed states. 

IV. ANTIFERROMAGNETIC INSTABILITY 
A. Reduction of Tn by local spin fluctuations 

According to eq. ()2.1H) . the Neel temperature Tn is determined by the competition be- 
tween the Kondo effect and the intersite exchange interaction: 

^ -^/s(0,Q)=0, (4.1) 



Xs(0) 4 

where Q is an ordering wave number to be determined. Effects of local spin fluctuations in 
SSA or those of dynamical mean fields in DMFT are included in the local term of l/x;s(0); 
/s(0, Q) corresponds to the conventional Weiss's mean field, which is an intersite effect and 
can only be included beyond SSA. 

The local susceptibility Xs(0) is almost constant at T ^ Tk, while it obeys the Curie law 
at T ^ Tk. We use an interpolation between the two limits: 

71 

Xs(0) = , (4.2) 



Here, Tk is an averaged one over Tk's at different sites. In disordered systems, the mapping 
conditions are different from site to site so that Tk's are also different from site to site. 
Such disorder in Tk causes energy-dependent lifetime widths, as is studied in Appendix. 
However, lifetime widths due to the disorder in Tk are small on the chemical potential 
in case of non-magnetic impurities. According to eqs. (|2.2j) and (j2.6|) together with the 
Fermi-liquid relation , ^^'^^ it follows that 

l/kBTK = 2W4p,{fx%^,, (4.3) 

in the absence of disorder. Then, a mean value of Tk in disordered systems is approximately 
given by eq. ()4.3|) with nonzero but small 7. It follows from eq. ()4.3|) that 

kBTK = ^-i-, (4.4) 

ctk 2Ws 

10 



with a numerical constant depending on n. As is shown in Fig. Efa), p-y{fi*) — 0.15 
for 0.1 ^ 7/1^*1 ~ 1 and 0.05 ^ \6\ ^ 0.25. We assume that is independent of n for 
the sake of simphcity: = 0.15. We are only interested in physical properties that never 
drastically change when c^k slightly changes. 



B. Exchange interaction arising from the exchange of pair excitations of quasi- 
particles 

The second term of eq. ()2.12j) is the sum of an exchange interaction arising from the virtual 
exchange of pair excitations of quasiparticles, Jqliuii, q), and the mode- mode coupling term, 
-4A(icU;,qr): 

2U^A7r,{iuJi, q) = Jc^{iuJu q) - 4A(i^/, q). (4.5) 
When higher-order terms in intersite effects are ignored, 



with 



jQiiUuq) = 4 



Xs(0) 



[Piiuji,q)-Po{iuJi)] 



P{iu,,q) = -j-J2g(^){ie^ + iu;,,k+q)g(^\ier.,k). 



(4.6) 



(4.7) 



e„ka 



The local contribution Po(i^/) = {^/^)J2qPi^^n,q) is subtracted because it is considered 
in SSA. The static component is simply given by 



P(0,g) = l^Z'r[e(fe + q)-^1-A[e(fc)-/i1 



e(fc)-^(fc + q)-iO 



(4.8) 



The magnitude of Jq(0, q) is proportional to fce^K or the quasiparticle bandwidth. 

Figures 0] and El show [P{0,q) — Po{0)] of the symmetric and asymmetric models. The 
polarization function at q = (±7r/a, ivr/a) is relatively larger in electron-doping cases than 
it is in hole-doping cases. 



C. Reduction of Tn by the mode-mode coupUng term 

Following previous papers,'^'^^'^^*^ we consider mode-mode coupling terms linear in in- 
tersite spin fluctuations Fs{iuji, q) given by eq. ()2.15|] : 

A(icui, q) = ki^i^iuJi) + Ki^oJu q) + A,(itU;, q). (4.9) 
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(a) (b) _ (c) 




q q q 



FIG. 4: Static polarization function [P{0,q) — Po(0)]|t*| of the symmetric model: (a) k-BT/\t*\ = 
7/[t*[ = 0.1, (b) kBT/\t*\ = 7/|t*| = 0.3, and (c) kBT/\t*\ = j/\t*\ = 1. Solid, dashed, broken, 
chain, and chain double-dashed lines show results for n ^ 1, n = 0.9, 0.8, 0.7, and 0.6, respectively. 
Here, F, X and M stand for (0,0), (vr/a, 0) and (vr/a, vr/a), respectively. 
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FIG. 5: [P(0,q) -Po(0)]|t*| of the asymmetric model: k^T l\t*\ = 7/|t*| = 0.1. A solid line shows 
a result for n ^ 1, dashed lines results of electron doping cases such as n = 1.1, 1.2 and 1.3, and 
broken lines results for hole doping cases such as n = 0.9, 0.8, and 0.7. 

The first term AL(iti;;) is a local mode-mode coupling term, which includes a single local 
four-point vertex function, as is shown in Fig. 4 of ref. Q. Both of K^iiuji, q) and Av(ico';, q) 
are intersite mode-mode coupling terms, which include a single intersite four-point vertex 
function; a single Fs{iuJiq) appears as the selfenergy correction to the single-particle Green 
function in As{iuJi,q) while it appears as a vertex correction to the polarization function 
in A^{iui,q), as are shown in Figs. 3(a) and 3(b), respectively, of ref. 
components are given by 

5 ksT 



4G. Their static 



Al(0) 



2Xs(0) N 



(4.10) 
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F^iiujuq) - ^J(qr)Xs(iwO 



(4.11) 



1 /c T 



2*.(0) iV ^^^^ 



with 



-[rf (k„)]%W(ie„ + io;o|, (4.13) 

5,(iu;i,p;g) = i^^^eT J]| ^ ^ (^("^(ie^, fc + g)^(°)(i£„, fc)^(°)(i^^ 

x^^f (ie„ + i.;,,fc + p) - [r(°)(ie„)]'[rW(ie„+i^O]'|, (4-14) 

with rCT°^(i£„) = {l/N)^f^g^a\i£n-,k). Because the selfenergy correction hnear in J[q) is 
considered in S IIIIl (1/4) J(q)Xs(ico'«) is subtracted in eq. ()4.1ip . 

In this paper, weak three dimensionahty in spin fluctuations is phenomenologically in- 
cluded. Because J{q) has its maximum value at g = (±7r/a, ±7r/a) and the nesting vector 
of the Fermi surface in two dimensions is close to q = (±7r/a, ivr/a) for almost half flUing, 
we assume that the ordering wave number in three dimensions is Q = (±7r/a, ±7r/a, iQz), 
with depending on interlayer exchange interactions. On the phase boundary between 
paramagnetic and antiferromagnetic phases, where eq. ()4.1|) is satisfled, the inverse of the 
susceptibility is expanded around Q and for small \uji\ in such a way that 

[1/Xs(iwi, Q + 9)]r=TN = ^(9) + + . (4.15) 

with 

A{q) = \A\^{q\iaf + [(g. - QM^ ■ (4-16) 

Here, c is the lattice constant along the z axis. Because Xs(i'^/5 Q + q) diverges in the limit 
of |q| ^ and cj^ on the phase boundary. 



5,(0, -Q;Q) = 5.(0, -Q;Q) 

-^b^eUE [9T(^^n,k-Q)f bi°)(ie„,fc)]'-[r(°)(i.„)]^ 1,(4.17) 



ml 

en (. k 
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can be approximately used for Bs{iuJi,p; q) in eq. ()4.1ip and B^{iui,p; q) in eq. ()4.12|) . Then, 
it follows that 



A(0,Q) 



2Xs(0) 



(4.18) 



with 



X 



16Vr^ 

^ " s 



de 



/;(e(fc)-/x*)+/;(e(fc+Q)-/x* 

(4.19) 
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In eq. (H?T^ . /^(e) is the derivative of f^{e) defined by eq. ^TT^ . In eq. (g^DI), P7^o(£^) is 
denoted by p{e), and 



P2(e) = Vp / de' 



J Pie') 



e-e' 



(4.22) 



In eq. 1)4.211) . the summation is restricted to |ti;| < ojc and |qr||| < gc- We assume that 
= n/3a and is given by a larger one of 8[t*\ and \J\. 

As was shown in the previous paper, ~ 1 for T/Ty^ ^ 1 and '-f/k-QTy: ^ 1. Physical 
properties for T/Tk ^ 1 or ■y/kBT}^ ^ 1 scarcely depends on a^. Then, we assume = 1 
for any T and 7 in this paper. 



D. Almost symmetric between 5 > and 5 < 

The Kondo temperature Tk is low in the limit of 5 ^ for nonzero T or 7. When 
T/Tk ^ 1 or 7/A;bTk ^ 1, the instability condition ()4.1|) becomes simple: 

1 5$ 1 

+ - 7^(Q) = 0- (4-23) 
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FIG. 6: Tn determined by eq. (|4.23() for J /\t\ = —0.3 and 5 = as a function of the exponent x 
of the anisotropy factor defined by |j42/j4||| = 10~^. 

When we ignore Al(0) = 5$/2xs(0) and we assume Tk = 0, eq. (j4.23|) gives Tn = 
(1/4) J(Q)/A;b = \J\/k-B- This Tn is nothing but Tn in the mean-field approximation for 
the Heisenberg model, to which the t-J model with the half filling is reduced. On the other 
hand, $ diverges at non zero temperatures for = 0; no magnetic instability occurs at 
nonzero temperature in complete two dimensions.— The reduction of Tn by Al(0) or by 
quasi-two dimensional thermal spin fluctuations depends on the anisotropy of Az/A\\. When 
only the superexchange interaction is considered, A\\ = \ J\. Although (1/4) Jq(0, Q + q) also 
contribute to the g-quadratic term, we assume A\\ = \J\ for the sake of simplicity. Figure IHl 
show Tn as a function of the exponent x of | A^/Ay | = 10^^ for the half filling (6 = 0). When 
we take x = 10, for example, we obtain Tn — 0.06|t|//cB or Tn — 0.2| J|//cb- We assume this 
anisotropy factor x = 10 in the following part of this paper. 

Figures [7| and |H1 show Tn of the symmetric and asymmetric models as a function of 
6 = 1 — n for various 7. The antiferromagnetic phase becomes narrower as 7 decrease. As 
long as 7 is almost symmetric between S > and 5 < 0, Tn is also almost symmetric between 
them. 

When 7 ^ W* or k^T ^ W*, Tn is determined by the competition between the superex- 
change interaction J and two quenching effects, l/xs(0) and Al(0). The reduction of Tn for 
small 1^1 is mainly due to Al(0). On the other hand, the critical concentration \6c\ below 
which antiferromagnetic ordering appears at T = K is determined by the competition 
between J and A;bTk, because thermal spin fluctuations vanish at T = K. Because the 
Kondo effect is rather weak, \6c\ is as large as \6c\ — 0.14, as is shown in Figs. [7| and |H1 
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6=l-n 

FIG. 7: Tn of the symmetric model as a fmiction of 5 = 1 — n. From the bottom, sohd, dashed, 
broken, chain, and chain double-dashed lines show Tn for 7/|t| = 0.2, 0.3, 0.4, 1, and 2, respectively. 
The topmost solid line shows Tn determined from eq. 1)4. 23(1 for comparison. The solid line with 
a positive slope shows ^/(p-y. Antiferromagnetic states whose Tn are much above and below this 
line are characterized as local- moment ones and itinerant-electron ones, respectively, according to 
discussion in SiVl 




FIG. 8: Tn of the asymmetric model as a function of 5. See also the caption of Fig. d Note that 
Tn is almost symmetric between 5 < and 5 > 0. 

When 7 < W* and k-QT < W* , well-defined quasiparticles can play roles in both of 
the enhancement and suppression of Tn. Antiferromagnetism is enhanced by jQ(i£:„,q), in 
which the nesting of the Fermi surface can play a significant role. However, the difference 
of the Fermi surfaces cannot give any significant asymmetry of Tn. The reduction of Tn 
by the Kondo effect is large. The intersite mode-mode coupling term, As(0, Q) + Av(0, Q), 
also suppresses Tn in addition to the local terms, /cb^k and Al(0). The quenching effects 
overcome the enhancement effect, and Tn decrease with decreasing 7. 
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V. DISCUSSION 



One of the most serious assumptions in this paper is that the homogeneous hfe-time 
width 7 of quasiparticles, which is phenomenologically introduced in this paper, can be 
defined. This is relevant when there are many weak impurities and their distribution is 
totally random. When disorder is so small that quasiparticles are well defined, electrons 
behave as Landau's Fermi liquid. The assumption is obviously relevant in such a case. On 
the other hand, when disorder is so large in the vicinity of the Mott transition or crossover 
that quasiparticles are never well-defined, electrons behave as localized magnetic moments. 
Because 7 plays no role, the theoretical framework of this paper must also be relevant in 
such a case. Then, an issues is whether the assumption is relevant or irrelevant in the 
crossover regime between the two regimes, for example, in the regime where the Anderson 
localization of quasiparticles occurs. Because the homogeneous lifetime width can be defined 
in the so called weakly localized regime of the Anderson localization,— we expect that the 
results of this paper apply to such a regime. The assumption is irrelevant when scattering 
potentials of dopants are strong and the concentration of dopants is rather small. The strong 
Anderson localization, where the localization length is short, is out of scope in this paper. 
The crossover between the Mott-Hubbard insulator and the Anderson insulator, which is 
examined in ref. ■32', is also out of scope in this paper. 

One may suspect that the theoretical framework of this paper cannot apply to the limit 
of (5 —s> or to vanishingly small \6\. According to the combination of Gutzwiller 's theory^ 
and the Fermi-liquid theory,'^*! the bandwidth of quasiparticles is vanishingly small but is 
still nonzero such as W* ~ 8\6t\. The Fermi liquid state is formed, at least, at T ^ 
8\6t\/kB', the spectral weight of quasiparticles is vanishingly small. It is obvious that when the 
renormalization of quasiparticles due to the Fock-type term of the superexchange interaction 
is taken into account the bandwidth becomes as large as \J\, at least, at T <^ 8|5t|//cB- 
We show in this paper that these quasiparticles are stable not only for k^T <^ 8\6t\ and 
7 ^ 8\6t\ but also for ksT ^ \J\ and 7 ^ \J\. The existence of elementary excitations 
or quasiparticles whose bandwidth is of the order of \J\ are implied by several studies on 
the t-J model on small lattices.^^SiMii^ Then, we can argue that when the Hilbert space is 
restricted within paramagnetic phases no divergence of the effective mass of quasiparticles 
or no M-I transition can occur except for J — 0, 5 ^ 0, and T K. What can occur for 
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nonzero J is a M-I crossover rather than a M-I transition, as is shown in Figs. and |21 

The coexistence phase of a metal and a Mott-Hubbard insulator is obtained for the 
Hubbard model with the just half filling and finite onsite U M It is interesting to examine 
whether or not the coexistence phase remains when the renormalization of quasiparticles by 
an intersite exchange interaction such as the superexchange interaction is included. It is also 
interesting to examine whether or not the coexistence phase remains for non-half fillings. 
Quasiparticles are well-defined for T ^ Tk and 7 ^ ^b^k- According to refs. I2JI and 



25l . the exchange interaction jQ(+iO,g), which arises from the virtual exchange of pair 



excitations of quasiparticles, shows an almost T-linear dependence consistent with the Curie- 
Weiss law at T ^ Tk in a small region of g ^ Q, with Q being the nesting wavenumber, for 
itinerant-electron antiferromagnets where the nesting of the Fermi surface is significant.^^ On 
the other hand, the T-linear dependence of l/xs(0) at T ^ Tk is responsible for the Curie- 
Weiss law of local-moment magnets or insulating magnets. Then, magnetism with Tn ^ 
Tk is characterized as local-moment one, while magnetism with Tn ^ Tk is characterized 
as itinerant-electron one. Local-moment magnetism appears for almost half fillings, and 
itinerant-electron magnetism appears for fillings substantially away from the half filling. 

Because no order is possible at non-zero temperatures in two dimensions^'^^ it is likely 
that Tn is substantially reduced by critical thermal fluctuations in quasi-two dimensions. 
As is examined in § HVDj the reduction of Tn is large; Tn — 0.2|J|/A;b for the anisotropy 
of lAz/Ayl = 10^^°. This explains observed Tn — 300 K for cuprate oxides, when we take 
\J\ ~ 0.15 eV. The exponent 10 seems to be little larger than actual ones. We should 
consider the reduction of Tn by thermal critical fluctuations more properly than we do in 
this paper. 

The nesting of the Fermi surface becomes less sharp as 7 becomes larger, as is shown 
in Figs. E] and El One may argue that antiferromagnetism is weaken by disorder in the 
itinerant-electron regime where the nesting of the Fermi surface plays a crucial role. How- 
ever, this effect is overcome by the reduction of Tk by nonzero 7. The renormalization of 
the quasiparticle bandwidth by the superexchange interaction becomes small with increasing 
7 so that the Kondo temperature /cbTk becomes small. On the other hand, the superex- 
change interaction is not reduced by nonzero 7. Then, antiferromagnetism is strengthen 
with increasing 7 even in the itinerant-electron region. 

An antiferromagnetic states appears in a narrow range of < \5\ ^ 0.02-0.05 in hole- 
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doped cuprates oxides {S > 0), while it appears in a wide range of < \S\ ^ 0.13-0.15 in 
electron-doped cuprate oxides (5 < 0). Tohyama and Maekawa^ argued that the asymmetry 
between hole-doped and electron-doped oxides must arises from the difference of the Fermi 
surfaces, and that the t-t'-J model should be used. They showed that the intensity of spin 
excitations is relatively stronger in electron doping cases than it is in hole doping cases. 
According to Fig. the polarization function at q = (±7r/2a, ±7r/2a) is relatively larger 
in electron doping cases than it is in hole doping cases. This asymmetry is consistent with 
that of spin excitations studied by Tohyama and Maekawa. However, the difference of the 
Fermi surfaces cannot explain the asymmetry of Tn, as is shown in Fig. |H1 

The condensation energy at T = K of the asymmetric model is also quite asymmetric;— 
it is consistent with the asymmetry of spin excitations and the polarization function dis- 
cussed above. On the other hand, Tn is significantly reduced by quasi-two dimensional spin 
fluctuations as well as local spin fluctuations of the Kondo effect. This large reduction of Tn 
arises from the renormalization of normal states; not only the Neel states but also param- 
agnetic states just above Tn are largely renormalized. It is plausible that the asymmetry of 
the condensation energy of paramagnetic states just above Tn is similar to that of the Neel 
states at T = K. It is interesting to examine by comparing the condensation energy of the 
Neel states and that of paramagnetic states just above Tn whether Tn is actually almost 
symmetric as is shown in this paper. 

Electrical resistivities of electron-doped cuprates are relatively larger than those of hole- 
doped cuprates are.-^ Then, we can argue that disorder must be relatively larger in electron- 
doped cuprates than it is in hole-doped cuprates, and that the asymmetry of Tn can arise, 
at least partly, from the difference of disorder between them. It is interesting to examine 
whether or not the symmetric behavior of Tn can be restored by preparing hole-doped and 
electron-doped cuprates with similar degree of disorder to each other. 

When nonmagnetic impurities of Zn ions are introduced on Cu02 planes, antiferromag- 
netic moments appear in the neighborhood of Zn ions.-^ Because the configuration of Cu ions 
is (3(i)^ and that of Zn ions is (3c?) the phase shift at the chemical potential due to scat- 
terings by Zn ions must be about tt/2 according to the Friedel sum rule. Then, scatterings 
by Zn ions must be very strong. One of the plausible scenarios is that the renormalization of 
/cbTk by the Fock-type term of the superexchange interaction is small in the neighborhood 
of Zn ions because of strong scatterings, that is, the quenching of magnetic moments by local 
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5=l-n 

FIG. 9: t* of the symmetric model as a function of 8 for various 7 and k-QT/\t\ = 0.02; Ws — 1 
or r = 0.5 is assumed instead of Ws — 2 or r = 1 (See text). From the top, sohd, dashed, broken, 
and chain hne show t* for ^ /\t\ = 0.04, 0.1, 0.2, and 0.4, respectively. For the sake of comparison, 
1/0^ is also shown by a bottom solid line. 




6=l-n 

FIG. 10: Tn of the symmetric model as a function of 5 for various 7; Ws ~ 1 or r = 0.5 is assumed 
instead of Ws — 2 or r = 1 (See text). From the bottom, solid, dashed, broken, and chain lines 
show Tn for 7/[t[ = 0.04, 0.1, 0.2, and 0.4, respectively. For comparison, Tn determined from 
eq. (|4.23() is shown by a topmost solid line. 

quantum spin fluctuations are weak there so that antiferromagnetism is enhanced there. It 
is interesting to develop a microscopic theory that can treat such an inhomogeneous effect. 

The three-point vertex function 0s and the mass-renormalization factor 0^ are those in 
SSA; Ws = 0s/07 — 2 as is discussed in § III Al When we take Ws — 2, on the other hand, 
the superexchange interaction constants as small as \J\ = 0.02-0.03 eV, which are much 
smaller than experimental ones \J\ = 0.10-0.15 eV, give as high as = 50-150 K ,^°i^^ 
which are as high as observed T^. This implies that the renormalization of the vertex 
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function and the mass renormalization factor by intersite effects such as antiferromagnetic 
and superconducting fluctuations should be included. Taking it into account, we treat Wg 
as another phenomenological parameter following the previous paper,— where Wg = 0.7-1 
is used instead of Ws ^ 2 in order to explain observed and observed coefficients of T- 
linear resistivities when \ J\ = 0.10-0.15 eV are used. Then, we replace Ws by rWg, with r a 
numerical constant smaller than unity; eq. ()3.5p is replaced by 



and eqs. ()4.3|) and ()4.4|) are replaced by those with rWg instead of H^s-~ Figures lUl and ITUl 
show t* and Tn, respectively, of the symmetric model as a function of 6 for r = 0.5 and 
various 7. The antiferromagnetic region extends with decreasing r. 

An antiferromagnetic state in the range of < 5 ^ 0.02 of hole-doped La2_5M5Cu04 (M= 
Sr or Ba) is characterized as a local-moment one. The so called spin-glass or Kumagai's 
phase^^ appears in the range of 0.02 ^ 6 ^ 0.05. We characterize the tail part of the solid 
line in Fig. ITUl as a critical line to the Kumagai phase;— 7/|t| — 0.04 are needed to reproduce 
Kumagai's phase. According to Fig. 121 \t*/t\ ~ 0.05 for 6 = 0.04. Then, we can argue that 
kpl ~ 2fcB7K/7 — 4|t*|/7, with kp the Fermi wave number and I the mean free path, must 
be 4-8 in Kumagai's phase. According to Fig. EJ^b), the nesting of the Fermi surface is 
substantial at least for 7/|t*| ^ 0.3; the nesting cannot be ignored for 7/|t*| ^ 1. Kumagai's 
phase must be a spin density wave (SDW) state in a disordered system rather than a spin 
glass. The divergence^ of the nuclear quadrupole relaxation (NQR) rate at Tn supports 
this characterization. 

VI. CONCLUSION 

The t-J model with J/\t\ = —0.3 on a quasi-two dimensional lattice is studied. First, an 
unperturbed state is constructed in SSA or DMFT. The Kondo temperature Tk or /cb^k is 
defined as the energy scale of local quantum spin fluctuations. The unperturbed state is a 
normal degenerate or almost degenerate Fermi liquid at T ^ Tk and is a non-degenerate 
Fermi liquid at T ^ Tk. No metal-insulator transition is possible within paramagnetic 
phases; what occurs is a metal-insulator crossover. The bandwidth W* of quasiparticles, 
which is about W* ~ 4A;bTk, is renormalized by the Fock-type selfenergy due to the su- 




(5.1) 
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perexchange interaction J, so that that W* is larger for lower T and smaller 7, with 7 being 
the phenomenologically introduced hfetime width of quasiparticles. The bandwidth W* is 
approximately given by the sum of terms of 0(|5t|) and 0(| J|), with 5 being the concentra- 
tion of dopants added to the just half-filled system; the term of 0(|5t|) is consistent with 
Gutzwiller's theory, and the coefficient of the term of 0(| J|) is 0(1) for /cBT/|iy*| ^ 1 and 
7/|H^*| ^ 1 and is vanishingly small for kBT/\W*\ > 1 or 7/|H^*| > 1. 

Next, antiferromagnetic instability of the unperturbed state is examined perturbatively 
in terms of intersite exchange interactions. This perturbation theory can treat both of 
local-moment magnetism and itinerant-electron magnetism. When Tn ^ Tk, with Tn the 
Neel temperature, or local thermal spin fiuctuations are dominant over local quantum spin 
fluctuations at Tn, magnetism is characterized as local-moment magnetism; it appear for 
almost half fillings of electrons. In the opposite case where Tn ^ Tk or local quantum ones 
are dominant over local thermal ones at Tn, magnetism is characterized as itinerant-electron 
magnetism; it appears for fillings away from the half filling. 

The reduction of Tn by quasi-two dimensional thermal critical fiuctuations is large when 
the anisotropy of the exchange interaction constants is large. For example, /cbTn/| J| 
in complete two dimensions, as is expected. The Neel temperature Tn — 300 K observed 
in cuprates with no doping can be explained, when we assume |J| = 0.10-0.15 eV and 
\Jz/ J\ — 10~^°, with J and the exchange interaction constants between nearest neighbors 
within a CUO2 plane and between CUO2 planes, respectively. 

When 7 > M^* or /cgT > W*, W* = 0{\St\) or /cbTr = 0{\5t\) so that the quenching 
of magnetic moments by local quantum spin fiuctuations is rather weak. Because thermal 
spin fluctuations vanish at T = K, an antiferromagnetic state is stabilized in a wide range 
of S. When 7 < W* and /cbT < W*, quasiparticles are well-defined. Because the nesting of 
the Fermi surface is substantial, the exchange interaction arising from the virtual exchange 
of pair excitations of quasiparticles is also responsible for antiferromagnetic instability in 
addition to the superexchange interaction. On the other hand, the quenching of magnetic 
moments by local quantum spin fiuctuations is strong. Not only quasi-two dimensional 
thermal critical fluctuations but also the quenching of magnetic moments by local quantum 
fluctuations make Tn substantially reduced or they destroy antiferromagnetic ordering. An 
antiferromagnetic state is stabilized only for small \S\; the critical concentration \Sc\ below 
which antiferromagnetic ordering appears is smaller for smaller 7. This result implies that 
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Tn is lower and a magnetic phase is narrower in a cleaner system. 

The difference or asymmetry of disorder between electron-doped and hole-doped cuprates 
must be, at least partly, responsible for that of Tn and antiferromagnetic phases between 
them. We characterize the so called Kumagai's phase as a spin density wave (SDW) state 
in disordered system rather than a spin glass. 
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APPENDIX A: SCATTERING POTENTIAL IN DISORDERED KONDO LAT- 
TICES 

In disordered Kondo lattices, the mapped Anderson models are different from site to 
site. The selfenergy for the jth one is expanded in such a way that T,j„{e + iO) = Sjcr(O) + 
^1 — e + - ■ ■ . When the energy dependence of the hybridization energy Aj(e) is ignored, 
it follows according to Shiba^ that 



with Edj and nja being the localized-electron level and the number of electrons with spin 
a, respectively. It follows from Eqs. ()3.(ij) and ()3.7|1 that cpj^ ^ (7r^/8)/|l — nj\ and (pjs — 
(77^/4)/! 1 — Ujl, with Uj = rij^ for almost half filling ^ 1. We assume non-magnetic 

impurities so that rij-^ = nj^. 

We denote the average number of electrons and their mean-square deviation by n = 
(nj)jjg and An^ = {{rij — respectively, where {■ ■ stands for an ensemble average 

over disordered systems. We assume that there is no correlation between different sites: 
{{rii - n){nj - n))^.^ = for i 7^ j. 

We consider the site-dependent part of E^j + + iO) as a static but energy-dependent 
scattering potential; it is approximately given by 



Edj + S,-.(+iO) = A,(0) tan [n (| 



(Al) 




(A2) 
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Here, the averages of 0j-y and Aj{0) are simply dented by 0-y and A. When we treat Vja{e) 
in the second-order SSA or the Born approximation, the coherent part of the ensemble 
averaged Green function is given by 

{G'^^\e + iO,k))^. =i , (A3) 

with ^ 

An'p,{e), (A4) 

for |5| < /cbTk- 

The bandwidth of quasiparticles is about 4kBTi<i; typical lifetime width is as large as 

ilmEf '^)(±A;bTk + iO) ~ -i{64/7T^)4>^An^\t\. (A5) 

07 

The energy-independent term can be ignored because 0^ ^ 1 for almost half fillings. It is 
quite likely that 0^An^ = 0(1) and '~f/\t\ = 0(1) for almost half fillings. 

It is straightforward to extend the above argument to a system in the presence of magnetic 
fields and a system with magnetic impurities. In such cases, riii — can be different from 
site to site; lmE''a^^\e + iO)/0-y can be large even for £ = 0. 



isf ^)(£+iO) = -i^ 
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